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Abstract 



- T— I I 

^ ■ We reconsider a cosmological evolution of domain walls produced 



by spontaneous breaking of an approximate discrete symmetry. We 
show, that domain walls may never collapse even if the standard bound 
on the vacuum energy asymmetry is satisfied. Instead of disappearing, 
these defects may form stable "bound states" - double wall systems. 
Possible stability of such a wall is a dynamical question and conse- 
quently restricts the allowed range of parameters. In particular, in 
the two Higgs doublet standard model with an anomalous Z2 symme- 
try, the above restriction suggests the mass of the pseudoscalar Higgs 
(would be axion) being close to the mass of the scalar one. 



1. Introduction 

Spontaneously broken discrete symmetries often play a fundamental role 
in particle physics models. Unfortunately, in the cosmological context such 
theories exhibit serious problems [1] as they lead to the formation of the 
topologically stable domain walls [2]. 

Enormous energy stored in these structures very soon comes to dominate 
the universe, unless the scale of discrete symmetry breaking (or the effective 
self coupling constant of the corresponding Higgs field) is very small (< 
lO'^GeV or so) [1,2]. 

One famous solution for the domain wall problem is inflation [4]. However, 
normally inflation deals with scales as high as say Mq ~ lO^^GeV (the GUT 
scale). In this situation, it is clear, that the topological defects formed below 
Mg, and in particular thous attributed to the electroweak phase transition, 
will survive. Therefore, it seems that at least for the weak scale domain walls 
the inflation may not be a good solution. 

Alternatively, cosmological troubles caused by the spontaneously broken 
discrete symmetry can be avoided if this symmetry is approximate, mean- 
ing that it is also explicitly broken by relatively small amount. From the 
naturalness point of view this 'solution' only makes sense if explicit break- 
ing results dynamically from some underlying physics e.g. from the anomaly 
[4]. Another possibility [5] is that gravity may not respect ungauged discrete 
(or continuous) symmetries and explicit breaking can manifest itself through 
the higher dimensional Planck scale (Mp) suppressed operators in the scalar 
potential. 



Whatever the source of the exphcit violation is, below certain temperature 
it creates energy difference between initialy degenerated vacua. Consequently 
one vacuum becomes energetically more favorable and corresponding domains 
start to expand pushing the walls away, whereas unfavorable domains tend 
to collapse and disappear. The necessary condition for the domain walls to 
disappear before dominating the universe has the following form [6] 

' > ^w/ « 

where e is the vacuum energy density difference and a is the wall energy 
per unit surface. In the simplest case, when the domain walls are formed by 
a single component real Higgs field $ with vacuum expectation value (VEV) 
< $ >= ±V", the condition (1) can be enough to guarantee that walls collapse 
without any trouble for cosmology. In the case of the real scalar field the 
expectation value < $ > always vanishes across the wall and two neighbor 
walls (bounding collapsing domain from the opposite sides) sooner or later 
will annihilate each other. However, in the realistic theories the Higgs fields 
usually are represented by complex scalars and have several components. 
The elements of the discrete symmetry group can be mimicked by certain 
abelian (or even nonabelian) phase transformation (which of course is not a 
part of continuous gauge group, since we are not interested in walls bounded 
by strings [7]). This circumstance brings a qualitatively new point in the wall 
structure implying that there is no topological reason for the absolute value 
of the Higgs field |$| to vanish inside the wall. In other words, whether this 
will happen becomes a matter of the parameter choice in the scalar potential. 



For a wide range of these parameters the energetically most favorable path 
in the field space, connecting two degenerated vacuums, never includes the 
point |<l>| = 0. Thus, the Higgs field never vanishes inside the wall and we 
can define its "winding number" through the defect. 

In the present paper we show that in such cases the domain walls may 
never collapse even if discrete symmetry is explicitly broken and (1) is sat- 
isfied. This has to do with the fact, that the walls of the collapsing domain 
can annihilate each other only if the Higgs field winds oppositely through 
this walls. Two walls with the same winding can never annihilate each other 
and disappear. Instead, they form classically stable "bound states"-double 
wall systems. Expectation value of the Higgs field returns to its initial value 
after traversing such composite wall. If winding is parametrised by phase 
arg$ = ^, we will have A^ = 27r through the wall. As far as we know, 
previously the 27r-walls have been considered only in the context of the walls 
bounded by cosmic strings [8], which often appear in the axion cosmology. In 
such case initialy there are no walls, but axionic string formed due to a global 
f/(l)pQ-symmetry breaking. 27r-walls are formed later on when U{1)pq gets 
explicitly broken (and only if there is no discrete subgroup left [9]). Each 
27r-wall is bounded by string. Existence of the string network crucially de- 
termines the cosmological evolution of the walls and makes them to decay 
without trace. 

In contrast 27r-walls discussed in the present paper can be cosmologically 
problematic, since there is no string network behind of their existence. They 
are formed as an bound states of would be topologically stable domain walls 
when laters collapse due to an explicit breaking of discrete symmetry. 



2. General mechanism 

In this section we will study mechanism of the 27r-wall formation on a 
simple model including one SU{3)— color triplet quark Ul,r transforming 
under discrete Z2 symmetry: 

ul -^ Ul, Ur -^ -Ur, (2) 

This quark gets a mass from the following Yukawa coupling 

Qu^ulUr + h.c (3) 

were $ = $/? + i$/ = |$|e*^ is a complex Higgs scalar ($ -^ — $ under 
^2), whose VEV breaks Z2 spontaneously. 

We choose a tree level scalar potential in the form: 

|$|2+_|$|4 ($2 + /i.c.) = |<|.|2 + _|$|4 |$|2cos2^ 



2 ' ' 4' ' 4 ' ' 2 ' ' 4' ' 2 

(4) 

For definiteness we assume all parameters m'^, M^, h being real and pos- 
itive. The possible quartic term $"^ which is allowed by all symmetries we 
have excluded for the simplicity of the analysis. Note that first two terms 
in (4) respect global f/(l)-symmetry $ -^ e*"$ which is explicitly broken by 
third one down to Z2 C U{1) and m is the mass of would be goldstone bo- 
son. The absolute value of the Higgs VEV is given by |$| = V = (M_±™_) 2 . 
This vacuum is discretely degenerated, since the potential is minimized by 
any 6 = ttN (where A^ is integer). This degeneracy results in the formation 
of the topologically stable domain walls. Topological constraint forces 6 to 



change by |A^| = n across the wall. However, this defects are not truly 
stable, because the Z2— symmetry is anomalous [4]. 

The important point about the wall structure is whether the absolute 
value of the scalar field vanishes inside the wall. Certainly, this is the case 
if m'^ > M^. In this situation the domain wall solution is given by kink 
(antikink) 

^R = ±Vtanh[xV{h/4:y/^]; (5) 

{x is a coordinate perpendicular to the wall) whereas the pseudoscalar com- 
ponent $/ is identically zero, since its effective [mass]^ term Mf = h^j^ + 
m'^ — M'^, which is spatially dependent trough the $r VEV, is positive every- 
where. In such a case, the wall network will evolve along the lines discussed 
in [4,5]. That is we expect that after Z2 will be explicitly broken by anomaly, 
walls can collapse without any cosmological harm (provided (1) is satisfied). 
The significant deviation from this scenario (to be discussed below) will 
occur if |$| stays nonzero everywhere including the vicinity of the domain 
wall. In this case the energetically most favorable trajectory connecting two 
neighbor vacuums 6 = txN and 9 = 71 {N ±1), goes through the saddle point 
$/ = ±[(M^ — 'm^)//^]-'^/^ and <I>k = 0, whereas the point |<I>| = is maximum. 
Of course, the condition M^ > m'^ is not yet sufficient to conclude that 
pseudoscalar will get a nonzero VEV inside the wall, since the gradient energy 
wants $/ to stay constant everywhere. So one has to check the stability of the 
small perturbation about $/ = in the background of kink. For example this 
can be done through the analysis of linearized equation for the pseudoscalar 
mode with the potential Mf on the existence of the bound state along the 



lines discussed by Witten [10] for the bosonically superconducting string. 
Such consideration shows that condensation of pseudoscalar inside the wall 
occurs at least for 

m^ <M^ ^ M^ - m^ (6) 

In such a case it costs a lot of energy (~ M^ per volume) for |$| to vanish 
somewhere in the space and in particular in the wall vicinity. To approximate 
the wall structure it is useful to put |$| = V = constant for a moment. The 
equation of motion for 9 becomes then a well known sine-Gordon equation: 

2||- sin 2^ = (7) 

(where ^ = mx and x is the coordinate transverse to the wall) which has a 
topologically stable solution (soliton) 

e{i) = 2-tan-^exp{i). (8) 

This is a reasonably good approximation for our purposes, although in 
reality |<l>| can not stay constant across the wall due to a back reaction. In 
turn this will alter the shape of d{^) , but whatever the exact form of 6{C,) 
is, it should change by |A^| = tt across the soliton. 

At the QCD scale the color instanton effects become important and they 
explicitly break Z2-symmetry. The effective Z2-noninvariant term generated 
in the scalar potential by instanton vertex has the form 

Vinst = -muA^ cos 6 (9) 



where A is a typical scale factor of the order the QCD scale Aqcd ~ 
200MeV and rriu is the quark mass. 

This term creates energy density difference e = 2muA.^ between 6 = 2N7r 
and 6 = (2A^+l)7r vacuums. This energy difference creates pressure difference 
which drives the domain walls. Assuming say rriuA^ > 0, the preference 
is given to ^ = 2N7r domains which tend to expand, whereas thous with 
6 = {2N + l)7r will tend to collapse and disappear. The dynamics of the 
collapsing walls is defined by several factors. At early stages (when explicitly 
violating effects are negligible) this are mainly wall tension and the frictional 
force. Crudely the time dependences of this two factors are proportional 
to ~ {aMp/t^y^"^ and ~ {Mp/ty respectively. At some point the system 
becomes dominated by pressure difference (with corresponding force per unit 
area ~ e) which forces the false vacuum to collapse. For walls to disappear, 
above has to happen before they dominate the universe and we are lead to 
the lower bound (1) for the vacuum energy asymmetry (for more details see 
[2,4,5,6]). 

Due to above arguments in general one expects [4] the domain walls dis- 
appearing soon after the QCD phase transition, since for the weak scale walls 
(destabilized by QCD anomaly) condition (1) is satisfied (e(Mp/o")^ ~ 10^° 
or so). 

The similar conclusion was made in [5] for the case of the domain walls 
destabilized by "gravity induced" higher dimensional operators. However as 
it will be shown in the present paper, the resulting dynamics of the walls 
is more sensitive to the parameters than one can naively expect from (1). 
Namely for the wide rang of parameters (forming the subrange of (1)) the 



walls may not collapse completely. In particular this can happen for (6) 
meaning, that the mass of the pseudoscalar mode (would be axion) is more 
than twice less the scalar one (radial mode). In such a case only the walls 
having opposite signs of A^ can annihilate each other, whereas others can 
not. Instead the laters will be paired forming stable double wall systems 
across which |A^| = 27r. Such a wall pairing can never occur if |$| vanishes 
across the wall. In this case the vacuas 6 = 2tiN (or 9 = {2N + l)7r ) 
corresponding to different A^ are identical. So in fact we have just two possible 
domains (say (+) and (— ) respectively). This has to do with the fact that 
actually 6 does not "winds" through the wall, since it becomes ill defined at 
|<I>| = 0. The planar domain wall solution separating (— ) and (+) ((+) and 
(— )) domains is given by kink (antikink) (5) rather then by soliton (8). Now 
if we switch on the instanton effects the (— ) domains will start to collapse and 
finally disappear, since the two neighbor domain walls (kink and antikink) 
will annihilate each other. 

Situation will drastically change if nonzero $ is energetically favored ev- 
erywhere in the space. That is, for example, if we are in the range of param- 
eters given in (6). Now 6 is well defined for any spatial point x (including 
thous in the vicinity of the domain wall) and as we show in the next section, 
this fact plays a crucial role in the subsequent evolution of the system. 

3. Cosmological evolution of |<I>| 7^ domain walls 

Now let us follow the history of the early universe in our model. At high 
temperature (T >> M, m) the ^-dependent part of the potential receives 
correction (Yukawa interaction is neglected) [11] 

9 



AV{T) = |$|'T'V6 (10) 

and we see that mass difference between the scalar and pseudoscalar fields 
(m^) never vanishes. Therefore, at the moment of the phase transition, when 
$R and $/ pick up a nonzero VEVs, they already "know" about explicit 
breaking of global would be f/(l)-symmetry. At this stage, probability of 
string-bounded wall formation is exponentially suppressed by the factor mR, 
where i? is a size of defect. Since the rang of our interest is m being not much 
(say within one order of magnitude) below M, the probability of large string- 
wall formation is negligible. So we are left with a network of domain walls 
of very complicated geometry. Traversing each wall, 6 winds by A^ = ivr 
continuously interpolating between two values that differ by tt. Difference 
among "winding numbers" of neighbor walls becomes very important after 
explicit breaking of the discrete symmetry. Instanton induced term (9) in the 
scalar potential gives preference to a certain domains (say to ^ = 2ttN) which 
tend to expand. Crucial point however is that not all 6 = (2A^+ l)7r domains 
will disappear. To see this consider collapsing domain 6 = n bounded by 
two infinite planar walls. It is clear that this boundaries can annihilate each 
other and disappear only if they correspond to the opposite windings of 6. 
That is two neighbor domains (of 9 = n one) must have both 6 = or 
both 9 = 27T (and never one and another 27r). Contrastly, if 7r-domain 
was initialy created between and 2tt ones, the corresponding walls can 
never annihilate each other since they have both the same A^. Instead of 
annihilation such walls will create stable bound states. Double wall solution 
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can be approximated analytically if we assume |$| = constant. In such case 
equation for 9 is double sine-Gordon equation. 



cPO 
2—- - sin2e - asinO = (11) 

where a = y2^^2 - The composite wall can be approximated by stable 
double-soliton solution 



e = 2arccos^^^—- — , Z^ = 2 + a/2 + 2,Jl + a/2th{2^ {I + a/2)i). (12) 

Note that the usual sine-Gordon equation (7) also admits 27r-soliton so- 
lution which is however unstable and tends to decay into two solitons with 
A^ = vr each. The last term in (11) oposses this decay and prevents two 
TT-solitons from escaping each other, forming the bound state. The width of 
the bound state (typical distance between A^ = vr walls) is 5 ^ —In-. Since 

\ J ^ 'ma 

A and V{^ m) correspond to QCD and weak scales respectively, this width is 
by order of magnitude larger the width of the constituent vr-solitons 5' ~ — . 
Infinite planar 27r-wall is not topologically stable and can decay quantum 
mechanically. This decay goes through the tunneling process which creates 
a hole connecting ^ = and 6 = 27i vacuums with each other. Imagine 
a closed path that starts in one of the domains (say in ^ = 0), goes once 
through the hole and finally returns back to the base point by traversing the 
wall in some other place. The one who travels along the path, will find phase 
6 changing by A^ = 27r at the end of the journey. Thus, |<l>| has to vanish at 
some point inside the region encircled by the path. By moving path around 
continuously (without crossing the edge of the hole), the point |<1>| = will 
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follow a closed line. So, there is a loop of the cosmic string enclosing the 
hole somewhere. Since, the space between two solitons is initialy filled with 
6 = 71 phase, the formation of the hole is costly in the energy of its own edge. 
The later includes the energy of the (cylindric) wall surface that bounds the 
hole and the energy of the string loop that is somewhere in this boundary. 
Hole will start to expand classically when, this energy will be overcomed by 
energy of the punched out piece of domain wall. The change in the energy 
caused by -R-radius hole formation can be estimated to be: 

AE = 7rR[2{jj^ + 6a)-R{6e + 2a)] (13) 

where a^ ~ 27rV'^ln— and a ^ AmV'^ are string and wall tensions re- 
spectively. Hole expands when its radius exceeds the critical value Re = 
{6 + fi'^/a)/2. The tunneling rate is exponentially suppressed by the factor 
2TTRla ~ TT[ln^{^f/^f{V/my. So that even for V/m -1-10 this wall 
can be stable for all practical purposes. Such 27r-wall of the horizon size (if 
formed) can be disastrous cosmologically. To destabilize this structures we 
should increase parameter rn^ and convert |$| = in to the saddle point. In 
the other words we have to increase the mass of angular pseudoscalar mode 
(would be goldstone boson) relative to scalar one (radial mode). Closer 
|$| = is to the saddle point, less costly is the string bounded hole forma- 
tion in the wall sheet. In the limit when |$| = becomes a saddle point, 
27r-walls become classically unstable. 

4. Two doublet extention 

Above model is trivially extendable to the realistic example of ref [4]. 

12 



This is the SU{2) U{1) standard model with two Higgs doublets $«, $d 
and anomalous discrete symmetry Z2 needed for natural flavor conservation. 
Under Z2: 

$« -^ -$«, $d -^ $d, Ql -^ Ql, Ur -^ -Ur, dR -^ dR (14) 

where Ql = ( , ) is the left-handed quark doublet. Z2-symmetry ensures 
that each Higgs doublet couples only to the fermions of the given charge. The 
scalar potential now takes the form 



V = E(^al*«r + ^^l-^^l") - ^l'^:'^^!' - h'\<^u?\^d? + Vphase (15) 

a 

Vphase = -\{K^df + h.C (16) 

where a = u,d. li h > the VEVs are pointing in the electrically neutral 
direction $a = ( " J . In this parameterization 



Vphase = -\Vyjcosi2e) 9 = 9^ - 9^ (17) 

For A > 0, Vphase is minimized hj 9 = nk {k = 1,2...) and there are topo- 
logically stable domain walls, li h + h' is larger than A, nonzero expectation 
values Vu,Vd 7^ will be energetically flavored everywhere in the space and 
even in the vicinity of the domain wall. So that as before, the winding of 9 
can be well deflned. Explicit violation of Z2-symmetry by QCD instantones 
leads to the formation of " composite" 27r- walls much in the way discussed in 
Sec. 2. The source of energy difference induced by instanton vertex in scalar 
potential has the form (for simplicity we assume just one family of quarks): 

13 



Vphase = -mufridA^cosO (18) 

where as before A ~ Aqcd — 200MeV and rriu, md are masses of -u and d- 
quarks respectively. The resulting double sine-Gordon equation for 9 has the 
form (11), where now a = "^yT^2 ■ Once again, the resulting 27r-walls can be 
practically stable and cause trouble for the cosmology, unless the parameter 
A is sufficiently large implying, that the mass of the would be axion is close 
to the Higgs mass. 

4. Conclusions 

We have shown that for the certain rang of parameters, compatible with 
standard bound, the collapse of unstable domain walls may end up with for- 
mation of classically (and practically) stable double wall systems. If formed, 
this structures can be problematic cosmologically and therefore their nonex- 
istence requires additional restriction of the parameter space. Of course, 
formation and subsequent evolution of this defects is very complicated pro- 
cess and requires much more detailed and perhaps numerical study. For 
example, collision of the domain walls may enhance the rate of hole forma- 
tion. At this point our analysis can be considered as 'ffist' approximation, 
but nevertheless, it provides an interesting enough example how cosmological 
considerations may restrict parameters in the particle physics models. 
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